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About the decomposition of pricing formulas 
under stochastic volatility models. 

Raul Merino*! Josep Vives! 


Abstract 

We obtain a decomposition of the call option price for a very general 
stochastic volatility diffusion model extending the decomposition obtained 
by E. Aids in [5] for the Heston model. We realize that a new term arises 
when the stock price does not follow an exponential model. The techniques 
used are non anticipative. In particular, we see also that equivalent results 
can be obtained using Functional ltd Calculus. Using the same general¬ 
izing ideas we also extend to non exponential models the alternative call 
option price decompostion formula obtained in [T] and [3] written in terms 
of the Malliavin derivative of the volatility process. Finally, we give a gen¬ 
eral expression for the derivative of the implied volatility under both, the 
anticipative and the non anticipative case. 


1 Introduction 

Stochastic Volatility models are a natural extension of the Black-Scholes 
model in order to manage the skew and the smile observed in real data. It is 
well known that in these models the average of future volatilities is a relevant 
quantity. Unfortunately adding a stochastic volatility structure, makes pricing 
and calibration more complicated, due to the fact that closed formulas not 
always exist and even when this formulas exist, in general, don’t allow a fast 
calibration of the parameters. 

During the last years different developments for finding approximations to 
the closed-form option pricing formulas have been published. Malliavin tech¬ 
niques are naturally used to solve this problem in [T] and [3] as the average future 
volatility is an anticipative quantity. Otherwise, a non anticipative method to 
obtain an approximation of the pricing formula is developed for the Heston 
model in [2]. The method is based on the use of the adapted projection of the 
average future volatility and obtain a decomposition of the call option price in 
terms of it. 
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In the present paper we generalize the results of [2] to general stochastic 
volatility diffusion models. Similarly, following the same kind of ideas, we extend 
the expansion based on Malliavin calcululs obtained in [T] and [3] . 

The main ideas developed in this paper are the following: 

• A generic call option price decomposition is found without the need to 
specify the volatility structure. 

• A new term emerges when the stock option prices does not follow an 
exponential model, as for example in the SABR case. 

• The Feynman-Kac formula is a key element in the decomposition. It 
allows to express the new terms that emerges under the new framework 
(i.e. stochastic volatility) as corrections of the Black-Scholes formula. 

• The decompostion found using Functional Ito calculus turns out to be the 
same as the decomposition obtained by our techniques. 

• We give a general expression of the derivative of the implied volatility, 
both for the non anticipative and the anticipative cases. 


2 Notation. 


Let S = {S{t),t S [0,T]} be a strictly positive price process under a market 
chosen risk neutral probability that follows the model: 



( 1 ) 


where W and B are independent Brownian motions, p G (—1,1), p: [0,T] x 
K+ —>■ M, 6: [0,T] x —>■ R+ and (j{t) is a positive square-integrable process 

adapted to the filtration of W. We assume on p and a sufficient conditions 
to ensure the existence and uniqueness of the solution of ©■ Notice that we 
don’t assume any concrete volatility structure. Thus, our decompositions can be 
adapted to many different models. In particular we cover the following models: 

• Black-Scholes model: p{t,S(t)) := rS'(f), 9{t, S{t),a{t)) := aS{t), p = 0, 
r > 0 and cr > 0. 

• CEV model: p{t,S{t)) := rS{t), 9{t, S{t),a{t)) := aS{t)^ with /3 G (0,1], 
p = 0, r > 0 and ct > 0. 

• Heston model: p{t,S{t)) := rS{t), 9{t, S{t),a{t)) := a{t)S{t), r > 0, 
a > 0 and 


df7^(t) = k(9 — (7^(t))dt + iyi/a^(t)dW(t), 


( 2 ) 


where k, 9 and v are positive constants satisfiyng the Feller condition 
2k9 > 
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• SABR model: fi{t,S{t)) := rS{t), 9{t,S{t),a{t)) := a{t)S{t)^ with /3 S 
(0,1], r > 0, CT > 0 and 

da{t) = aa{t)dW(t) (3) 


with a > 0. 

For existence and unicity of the solution in the Heston case see for example 
HU, Section 2.2. For the CE\^ and SABR models see and the references 
therein. 

The following notation will be used in all the paper: 

• We will denote by BS{t,S,a) the price of a plain vanilla european call 
option under the classical Black-Scholes model with constant volatility cr, 
current stock price S', time to maturity t = T — t, strike price K and 
interest rate r. In this case, 

BS {t, S, a) = S^{d+) - Ke-^^<^{d-), 

where $(•) denotes the cumulative probability function of the standard 
normal law and 

ln(S/R) + (r±4)r 

d± = - 1= -^—. 

<y\/T 

• We use in all the paper the notation Et[-] := where {Tt,t > 0} is 

the natural filtration of S. 

• In our setting, the call option price is given by 

R(t) = e-’’"Et[(S(r)-R)+]. 

• Recall that from the Feynman-Kac formula, the operator 

satisfies CgBS(t, S{t), 6{t, S(t), cr(t))) = 0. 

• We will also use the following definitions for j/ > 0: 

G(t,S(t),y) := S^{t)dlBS{t,S{t),y), 

Hit,S{t),y) ■.= S{t)dsG{t,Sit),y), 

K{t,S{t),y) := S\t)dlG{t,S{t),y) 

and 


L{,t,S{t),y) 


0{t,S{t),y) 

S{t) 
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3 A decomposition formula using Ito Calculus. 

In this section, following the ideas in [5], we extend the decomposition for¬ 
mula to a generic stochastic volatility diffusion process. We note that the new 
formula can be extended without the need to specify the underlying volatility 
process, obtaining a more flexible decomposition formula. When the stock price 
does not follow an exponential process a new term emerges. The formula proved 
in [2] is a particular case. 

It is well known that if the stochastic volatility process is independent from 
the price process, the pricing formula of a plain vanilla European call is given 

by 


V{t)=Et[BS{t,S{t),am 

where (t) is the so called average future variance and it is defined by 

Naturally, a{t) is called the average future volatility. See [TU], pag. 51. 

The idea used in [5] consists in consider the adapted projection of the average 
future variance 


:= Et(cr^(t)) = ^7^ Et[cr2(s)]fis. 

and obtain a decomposition of V (t) in terms of v(t). This idea switches an antici- 
pative problem related with the anticipative process a{t) into a non-anticipative 
one with the adapted process v{t). We apply this technique to our generic 
stochastic differential equation (HD. 

Theorem 3.1. (Decomposition formula) For all t G [0,T) we have 


Vit) = BS{t,S{t),v{t)) 


-Et 

2 


rE* 


2 


e *'^G{u,S{u),v{u)) (^L^{u,S{u),a{u)) — du 


e *^K{u,S{u),v{u))d[M,M]{u) 


e *^'^L{u,S{u),a{u))H{u,S{u),v{u))d[W,M]{u) 


where M(t) := Jq Et [o'^(s)] ds = f* a^{s)ds + {T — 

Proof. Notice that BS{T, S{T),v{T)) = e~‘^'^V{T). As e“’’W(f) is a mar¬ 
tingale we can write 

e-"V(t) = Et {e-^^V{T)) = E* {e-^^BS{T, S{T),v{T))) . 

Our idea is to apply the Ito formula to the process e~^^BS{t, S{t), v{t)). 
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As a consequence of the fact that the derivatives oi BS are not bounded we 
have to use an aproximation to the identity argument changing BS(t, S, a) by 

BSnit,S,a) := a)V-n(5) 

where ipn{S) = for some (j) € such that ^(5) = 1 for all IS”! < 1 and 

0(5') = 0 for all |5| > 2, and v{t) by (^S + E [(T2(s)ds]^ , where 

e > 0, and apply finally the dominated convergence theorem. For simplicity we 
skip this molliher argument in all the paper. 


So, applying the Ito formula, using the fact that 

daBS{t, 5, a) = S^ardjBSit, 5, a) 

and the Feynman-Kac operator (|3]), we deduce 

e-^^BS{T, 5(T), v{T)) - e-’'‘B5(t, S{t), v{t)) = 

= [ e-^'^£ysBS{u,S{u),v{uj)du 


( 5 ) 


e '~'^dsBS{u, S{u),v{u))9{u, S{u),a{u)) (^pdW{u) + \/l — p^dB{u)^ 


+ 


e-™5^(u)a|B5(M, S{u),v{u))dM{u) 


+ 




e '''^S^{u)dgBS{u, S{u),v{u)) [L'^{u, S{u),a{u))du — a‘^{u)du\ 
e"™ {S^{u)dl {S^{u)dlBS{u, S{u), v{u)))) d [M, M] {u) 

S{u), a{u)) (9s (5 ^(u)9|B5(u, 5(m), v{u)))) d [W, M] (u). 


Taking conditional expectation and multiplying by e'"*, we have: 

Etle-’'(^-*^BS(T, S(T), v(T))] = BS(t, S(t),v(t)) 

1 , 




e *^S^{u)dgBS{u,S{u),v{u))[L^{u,S{u),a{u))du — a^{u))du 


rEt 


-Et 

2 


o-r{u-t) ( c-2 


(5^(m)9| {S‘^{u)dlBS{u, 5(w), u(w)))) d [M, M] (u) 


e *'’0{u,S{u),a{u))ds {S‘^(u)d‘gBS{u,S{u),v{u))) d\W,M](u) 


□ 


Remark 3.2. In the following operators are defined for X{t) — log5(t) 
. G{t,X{t),a{t)) := {dl-d,)BS{fiX{t),a{t)). 

. H{t,Xit),a{t)) := {dl-dl)BS{fiX{t),a{t)). 
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. k{t, X{t),a{t)) := (54 _ 293 + dl) BS{t, X{t),a{t)). 

We observe that 

• G{t,X{t),a{t)) = G{t,S{t),a{t)). 

• k{t,X{t),a{t)) = K{t,S{t),a{t)). 

• H{t,X{t),a{t)) = H{t,S(t),a{t)). 

Remark 3.3. We have extended the decomposition formula in m to the generic 
SDE ([T|) . When we apply the ltd calculus, we realize that Feynman-Kac formula 
absorbs some of the terms that emerges. Finally, we ended up with three new 
terms to adjust the price. It is important to note that this technique works for 
any payoff or any diffusion model satisfying Feynman-Kac formula. 

Remark 3.4. Note that when 6{t, S{T),a{t)) = a{t)S{t) (i.e. the stock price 
follows an exponential process) then 


V{t) = BS{t,S{t),v{t)) 


+ 77^4 




S{u), v{u))d [M, M] (w) 




a(u)H(u, S{u), v(u))d [W, M] (u) 


and the term 

rT 


-Et 

2 


e S^{u)dgBS{u, S{u),v{u)) [L^{u, S{u),a{u)) — du 


vanishes. 

Indeed, we will show that due to the use of Feynman-Kac formula this is 
happening. 


Movement of the asset + Movement of the volatility 
= Su)'^dgBS{u, S{u), v{u))du + dcBS{u, S{u), v{u))dv{u) 

= ^a'^{u)S^{u)dgBS{u, S{u),v{u))du 

+ S{u),v{u)){dM + v'^du — a'^du) 

= ^S'^{u)dgBS{u, S{u),v{u)){dM + v'^du) 

where ^S^{u)dgBS{u, S{u),v{u))v^ is used into the Feynman-Kac formula and 


E. 


-S^dlBS{u, S{u),v{u))dM 


= 0 . 
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4 Basic elements of Functional Ito calculus. 

In this section we give the insights of the Functional Ito calculus developed 

in [3E1I71IH1- 

Let X : [0, T] X i — > K be an Ito process, i.e. a continuous semimartingale 
defined on a filtered probability space (17, P) which admits the 

stochastic integral representation 

^(u)du + f a(u)dW(u) ( 6 ) 

Jo 

where W is a Brownian motion and fi{t) and a(t) are continuous processes 
respectively in x [0,T]) and L^{il x [0,T]). 

We define i7([0,T],K) the space of cadlag functions. Given a path x G 
i7([0,T],R), we will denote as x* its restriction to [0,f]. For h > 0, the hori¬ 
zontal extension Xt^h is defined as 


Xt.hiu) = Xtiu) = xiu), uG[0,t[; 

and the vertical extension as 

Xt.hiu) = xit), uGit,tYh] 

(7) 

Xt('u) = Xtiu) = xiu), 

u G [0,4; 

(8) 


Xtit) = x{t) + h, i.e. Xt{u) = x{u) + hl^t=u}■ 

A process Y : [0, T] x K —>■ K, progressively measurable with respect the natural 
filtration of X, may be represented as 

Yit) = F {t, {A(s), 0 <s<t})= Fit, Xt) 

for a certain measurable functional F : [0,<:] x i7([0,t],R) —>■ R. Let F°“ be the 
space of locally lipschitz functionals with respect the norm of the supremum on 
i7([0, f + /i], R), that is, it exists a constant C > 0 such that for any compact K 
and for any x G 77([0, t],K) and y G 77([0, t + h],K) we have 


X(t) = xo + 


17^(7, Xt ) F it Y h, Y 11 oo ■ 

Under this framework, we have the next definitions of derivative: 


Definition 4.1. (Horizontal Derivative) The horizontal derivative of a func¬ 
tional F G F°“ at t is defined as 


VtFit,Xt) 


Fit + h,xt,h) - Fit,xt) 

lim -;- 

/t-s-o+ h 


(9) 


Definition 4.2. (Vertical derivative) The vertical derivative of a functional 
F G F°“ at t is defined as 


VxFit,xt) 


lim 


Fit,x’() - Fit,xt) 
h 


( 10 ) 


Of course we can consider iterated derivatives as V^x- 

We also have the following ltd formula that works for non-anticipative func¬ 
tionals: 
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Theorem 4.3. (Functional ltd Formula) For any non-anticipative functional 
F S F°“ and any t G [0,T] we have 



'DuF{u,Xu)du+ / y xF{u,Xu)dX{u) 


provided T>tF, X^F and XxxF belong to F°“. 
Proof. See [HlIZ]. 


□ 


5 A general decompostion using Functional Ito 
Calculus. 

In this section, we apply the technique of functional Ito calculus to the 
problem of finding a decomposition for the call option price. The decomposition 
problem is an anticipative path-dependent problem, using a smart choice of the 
volatility process into the Black-Scholes formula we can convert it into a non 
anticipative one. It is natural to wonder if the functional ltd calculus brings 
some new insides into the problem. 

We consider the functional 


F{t, S{t),a() = e-^*BS{t, S{t),f{t, a^)) 


where is the path-dependent process and / G F°“ is a non-anticipative func¬ 
tional. 

Under this framework, we calculate the derivatives using the functional Ito 
calculus respect the variance and we write them in terms of the classical Black- 
Scholes derivatives. We must realize that, for simplicity, the new derivatives are 
calculated respect the variance instead of the volatility of the process. 

Remark 5.1. If d denotes the classical derivative, we have: 

• Alternative Vega : 


X,.F = e-^*dfBS{t, S{t), fit, a()) X,.f{t, a^). 


• Alternative Vanna: 


'^a^.sF = df^sBSit,S{t),fit,a'()) V<,2/(t,cr|). 


• Alternative Vomma: 


= e-'-^dfjBSit,Sit),fit,a()){X,.fit,a()f 
- e-^^dfBSit, S{t), fit, a())Xl.fit, a(). 


• Alternative Theta: 


VtF = -re--*BSit,Sit),fit,a()) 

+ e-^^dtBSit,Sit),fit,a()) 

+ e~^*dfBSit,Sit),fit, a1))Vtfit, a")). 


Theorem 5.2. (Decomposition formula) For all t G [0,T), S{t) and f{t, a^) > 0 
we have 


V{t) = BSit,S{t)Jiu,a()) 


+ Et 

+ -Et 
2 

+ -Et 
2 

+ pEt 


al)TG{u,S{u), f{u, al))Vuf{u, al)du 
e“’'(““‘)G(w, 5'('u), /(w, al)) {L^{u,S{u),au) - f{u, du 
al fT^K{u,S{u),f{u, al))d [f{u, al)] 

r 

5'(w), a{u))H{u,S{u), f(u, (jl))f{u, aDrd [Wal)] 


Proof Notice that F{T,X{T),a^) = e-^^BS{T, SiT), f{T,a^)) = e-'^’^Vr- 
As e~'^*'V{t) is a martingale we can write 

e-’-V(t) = Et {e-^^V{T)) 

= Et{e--^BS{T,SiT),f{T,a^r))) 

= Et(F(T,5(T),4)). 

Our idea is to apply aproximation to the identity argument as in Theorem 13.11 
and then use the functional Ito formula to 


F{t, S{t), 4 ) = e-^*BS{t, S, fit, 4 )). 


We deduce that 

FiT,SiT),al)-Fit,Sit),a() 

= 'DuFiu,Siu),al)du + J^ VsFiu, Siu),al)dSiu) 

+ V,.Fiu, S{u), al)du + i V|F(t, S{t),<7^^)d[S, 4(u) 

+ ^ '^l^Piu,Siu),al)d[al,al] VIF{u, Siu),(jl)du 

+ j\l,.Fiu,Siu),al)d[S{u),<yl] 

+ ^ df{"^sFiu,Siu),al))d[Siu),f{u,al)\ 

+ ^ df {"^a-^Fiu,S{u),al))d[a‘^,fiu,al)\. 


Note that: 

• As S{t) is not path-depeirdent, we have that V 5 (-) = ds{‘)- 

• As u > t and / is a noir-airticipative functional, then Vcr 2 („)/(t, 4) = 0- 
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So, we have 


F{T,S{T),a^)-F{t,Sit),af) 


'DuF{u,S{u),al)du + dsF{u, S{u),al)dS{u) 
dlF{u,S{u),al)d[S,S]{u) VlF{u, S{u),al)du 


dj,sF{u,S{u), al)d [S{u),f{u, al)] . 


We deduce that 


F(r,5(T),4)-F(t,5(t),a,^) 

= ^ Cf^^^,2jBSdu +e-'^^dfBS{u,S{u),f{u,al))V^f{u,al)du 

+ ^ e~'^^d‘^BS{u,S{u),f{u,al)){9'^{u,S{u),a{u))-S'^f{u,al))du 

+ dsBS{u,S{u),f{u,<7‘^))9{u,S{u),<7{u)){pdW{u)F\/^^^^^dB{u)^ 

+ \ e-^''djBS{u,S{u),f{u,al))d[f{u,al),f{u,al)] 

+ P e~^''d‘jsBS{u,S{u)J{u,al))9{u,S{u),a{u))d[W{u)J{u,al)\. 

Taking now conditional expectations, using and multiplying by e’’* we obtain 
that 


e-(T-t)E, [F{T,S{T),al)] = BS{t, S{t), f{u,af)) 




£ e-™/(u, 4)(T - t)S^dlBS{u, S{u), f{u, 4))P„/(u, 4) 


du 


-Et 

2 


+ -Et 
2 


+ pEt 


e ™G(w, S'(u), /(w, 4)) {L^{u,S{u), a(u)) - f(u, 4)) du 


e ™/(m, crl)^T^K{u,S{u), f{u, al))d [f{u, 4), /(«, 4)] 


e ™i?(w, S'(u), /(u, crl))f{u, al)T9{u, S{u),a{u))d [W{u),f{u, 4)] 


□ 

Remark 5.3. Note that functional ltd formula proved in m holds for semi¬ 
martingales, but in m is also proved for Dirichlet process. In both cases, the 
hypothesis hold by definition of f and differentiability of the derivatives of Black- 
Seholes function when T,S,a > 0. Therefore, this technique can be applied to 
these models. 
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Remark 5.4. Note that Theorem \5.S\ coincides with Theorem no when we 
choose the volatility function as f(t,a^) = v(t). We found an equivalence of 
the ideas developed by and m in the decomposition problem. Both 

formulas come from very different places, the ideas under are based 

on an extension to functionals of the work while the main idea of m is to 
change a process by his expectation. Realise that standard Rd calculus also can 
be applied to Dirichlet process (for more information see f3^). 

Remark 5.5. Realize that Theorem \5.‘A holds for any non-anticipative f(t,a(). 
R is no trivial to find a different non anticipative process f(t, uf) different from 
the one chosen in m- 

6 Basic elements of Malliavin Calculus. 

In the next section, we present a brief introduction to the basic facts of 
Malliavin calculus. For more information, see m- 

Let us consider a Brownian motion W = {W{t),tG [0,T]} defined on a 
complete probability space (fl, P). Set H = L^([0,T]), and denote by W{h) 
the Wiener integral of a function h G H. Let S be the set of random variables of 
the form F = f{W{hi,... ,W{hn)), where u > 1, / G C“, and hi,... ,hn € H. 
Given a random variable F of this form, we define its derivative as the stochastic 
process F,t G [0,T]} given by 

n 

Df^F = J2d.JiW{hi),...,W{h,,))h,{x), te[0,T]. (11) 

The operator and the iterated operators are closable and unbounded 
from Lf{Tt) into L^([0,T]"’ x 11), for all n > 1. We denote the closure of S with 
respect to the norm 

n 

:= Wlh,!,, + E ■ (12) 

k=l 

We denote by the adjoint of the derivative operator . Notice that 
is an extension of the Ito integral in the sense that the set L^([0,T] x H) of 
square integrable and adapted processes is included in Dom6 and the operator 
5 restricted to L^([0,T] x H) coincides with the Ito stochastic integral. We use 

the notation S{u) = u{t)dW{t). We recall that := Lf is 

contained in the domain of 5 for all u > 1. 

We will use the next Ito formula for anticipative processes. 

Proposition 6.1. Let us consider the processes X{t) = x{0) + u(s)dW(s) + 
fg v(s)ds, where u,v G L^([0,T] x ^l). Furthermore consider also a process 
Y(t) = 9{s)ds, for some 9 G L^’^.Lel F : ^ R a twice continuously 

differentiable function such that there exists a positive constant C such that, for 
all t G [0,T], F and its derivatives evaluated in (t, X{t),Y{t)) are bounded by 
C. Then it follows 

F{t,X{t),Y(t)) = F(0,X(0),F(0))+ / dsF{s,X{s),Y{s))ds 

Jo 
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+ [ d^F{s,Xis),Y{s))dXis) 

Jo 

+ [ dyFis,Xis),Y{s))dYis) (13) 

Jo 

+ [ dl^yF{s,X{s),Y{s)){D-Y){s)u{s)ds 

Jo 

+ II dlF{s,Xis),Y{s))u\s)ds, 
where {D~Y){s) := D^Y(r)dr. 

Proof. See [T]. □ 

The next proposition is useful when we want to calculate the Malliavin 
derivative. 


Proposition 6.2. Let a and b be continuously differential functions on R with 
bounded derivatives. Consider the solution X = {Xt,t G [0,T]} of the stochastic 
differential equation 

(T(X(s))(ilT(s) + f b{X{s))ds. 

Jo 

Then, we have 


X(t) = x{0) + 


DsX{t) 


a(X(s)) exp 



a'(X(s))dW(s) 



where A(s) = [b' — ^ (cr')^](Ar(s)). 

Proof. See [T^], Section 2.2. □ 


7 Decomposition formula using Malliavin calcu¬ 
lus. 

In this section, we use the Malliavin calculus to extend the call option price 
decomposition in an anticipative framework. This time, the decomposition for¬ 
mula has one term less than in the Ito formula’s setup. 


We recall the definition of the future average volatility as 
d(t) := J — - f a‘^{s)ds. 


T-t.. 

Theorem 7.1. (Decomposition formula) For all t G [0,T), we have 
V{t) = V,t[BSlf,S{t),d{t))] 


2 


e ''‘^G{u,S{u), cr(u)) {Lf(u, S{u), a{u)) — cr^(u)) du 
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where 


rE* 


J e *'^L{u,S{u),a(u))H(u,S{u),au) D'^a^(r)d^du 


and 


Git,S{t),a{t)) := S^{t)diBS{t,Sit),cr{t)) 
H{t,S{t),a{t)) := Sit)dsG{t,S{t),a{t)), 


L(t,S(t),a(t)) := 


S{t) 


Proof. Notice that e BS{T, S{T),a(T)) = e As e ^*V{t) is a martin¬ 

gale we can write 

e-"V(t) = Et {e-^^V{T)) = E* {e-^^BS{T, S{T),d{T))) . 

So, using the aproximation to the identity argument and then applying the Ito 
formula presented in Proposition 16.11 to 

e-^*BS{t,S{t),d{t)). 

We deduce and using ([S]) and dl]) that 

e-^^BS{T, S{T),d{T)) - S{t), d{t)) = 

_ f CgsBS{u, S{u),d{u))du 


1 


e ^'^S^{u)dgBS{u,S{u),(j{u)) 


0(u, S(u), a(u)) 


S(u) 


du 


e ^'^S^{u)dgBS{u, S{u),a{u))a^{u)du 


dsBS{u, S{u),au)d{u, S{u), (j(u)) (^pdW(u) + \/l — p^dBiu^ 

+ ^ e~^'^e{u,S{u),a{u))ds {S'^{u)djBS{u,Siu),d{u))) (j D'^a'^{r)dr 


+ I ^ 

pT 


Taking conditional expectation and multiplying by we have 
Et[e-’'(^-‘)B5(r, S{T), d{T))] = E* [BSit, S{t), a(t))] 


+ -Et 
2 


+ fE. 


e *''G{u,S{u),a{u)) {^L'^{u,S{u),a{u)) — cr'^{u)) du 

L{u, S{u),a-{u))H{u, S{u),d{u)) ( J (r)dr'\ du 


□ 


du. 
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Remark 7.2. As it is expected, a new term emerges when it is considered ([T]) 
like it happen in Theorem \3.1[ 

Remark 7.3. In particular, when 9{t, S{T),a{t)) = a{t)S{t) 


V{t) = BS{t,StMt)) 


+ fE. 


e ^^a{u)H{u,S{u),au) ij D^a^{r)drjdu 


Also, the gamma effect is cancelled as we have seen in the ltd formula section. 

Remark 7.4. Note that when v(t) is a deterministic function, we have that all 
decomposition formulas are equal. 

Remark 7.5. When p = 0, we have 

Et [BS{t, Sit), ait)) - BSit, St, ^;(^))] 


= -Et 
2 


- -Et 
2 


- -Et 






(G'(u, Siu), cr(u)) — Giu, Siu),viu))) L^iu, S'(u), aiu))du 


(G'(u, Siu), cr(u)) — Giu, Siu),viu))) a'^iu)du 


e *'^Kiu,Siu),viu))d[M,M]iu) 


In particular, when 9 it, Sit), ait)) = a it) Sit) : 

Et [BSit, Sit), ait)) - BSit, St,vit))] 


= -^Et 


e-A^-*)Kiu, Siu), viu))d [M, M] (w) 


The difference between the two appoaches is given by the vol vol of the option. 


8 An expression for the derivative of the implied 
volatility. 

In this section, we give a general expression for the derivative of the implied 
volatility under the framework of Ito calculus and Malliavin calculus. There 
exist a previous calculation of this derivative in the case of exponential models 
using Malliavin calculus in [3]. 


Let /(S'(t)) denote the implied volatility process, which satisfies by definition 
Vit) = BSit, Sit), I is it))). We calculate the derivative of the implied volatility 
in the standard Ito case. 


Proposition 8.1. Under for every fixed t G [0, T) and assuming that 
(u(t))“^ < c» a.s., we have that 


dsIiS*it)) 


Et 


dsF 2 iu, S*iu), viu))du 


d,,BSit,S*it),IiS*it))) 
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where 


E, 


{Fi{u,S*{u),v{u)) + dsFs^u, S*{u),v{u))) du 


2Sd,BSit,S*it),I{S*{t))) 


E, 


= -Et 
2 


Fi{u, S{u), v{u))du 


e *'^G{u, S{u),v{u)) (^L'^{u, S{u),a{u)) — (j'^{u)) du 


rE. 


+ -Et 
2 


S{u),viu))d [M, M] (u) 


,-r{u- 


t) 9{u,S{u) a{u)) 5(u), z;(u))d [W, M] {u) 


S{u 


and 


Et 


= -Et 
2 


Et 


= -Et 
2 


rEt 


F 2 {u, S{u),v{u))du 


^-r{u- 


t) e{u,s{u) Mu)) 


S{u) 


F^^u, S{u), v{u))du 


e *M{u, S{u),v{u)) (^L'^{u, S{u),a{u)) — a'^M) du 


^-r(u-t) 


K{u, S{u),v{u))d [M, M] (u) 


Proof. Taking partial derivatives with respect to S{t) on the expression V{t) = 
BS{t, S{t), I(S{T))), we obtain 

dsV{t) = dsBS{t, Sit), liSiT))) + d^BSit, 8(1), IiSiT)))dsI{S{t)). (14) 
On the other hand, from Theorem 13.II we deduce that 


Vit) = BSit,Sit),vit))+Et 
which implies that 

dsVit) = dsBSit,S{t),vit)) + Et 


Fiiu, S{u), v{u))du 


dsFiiu, S{u), v{u))du 


(15) 


(16) 


Using that (u(t)) ^ < oo we can check that dsVit) is well define and hnite a.s. 
Thus, using that S'*(t) = K exp(r(T — t)), (fHl) and (fTBl) . we obtain 


dsiMit)) = 


dsBSjt, S*it), vjt)) - dsBSjt, S*it), IjSjt))) 
d,BSit,S*it),IiSit))) 
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ft dsFi(u,S*(u),v(u))du 


d^BS(t,s*(t),i(s(t))) • 

From m we know that = 0, where /°(t) is the implied volatility in the 

case p = 0, so 


dsBS{t,S*it),v{t)) = dsBS{t,S*{t),I°iSit)))-Et 
So, we have that 


dsFsiu, S*{u),v{u))du 


dsI{S*{t)) = 


dsBS{t, - dsBS{t, 


E, 


d^BSit,S*{t),IiS*{t))) 


ft dsF 2 {u,S*{u),v{u))du 


d^BSit,S*it),I{S*it))) 
On the other hand, we have that 

dsBS{t,S*{t),v{t))=m 

and 


BSit,S*it),vit)) = S im - 

where (j) is the standard Gaussian density. Then 

BS{t,S*{t),v{t)) + S 


dsBSit,S*{t),vit)) = 


2S 


and 


dsBS{t, S*it),I°it)) - dsBS{t, S*(t),I{S*(t))) 
^ {BS{t, - BS{t, 


= -Et 

2S 


{Fi{u, S*{u),v{u)) + dsF 3 {u, S*{u),v{u))) du 


□ 


Now, we derive the implied volatility using the Malliavin calculus. It has 
been proved in [2, in the case when 6{t, S(T), a{t)) = a{t)S{t). 

Proposition 8.2. Under o, for every fixed t G [0,r), assuming that (ft(t)) ^ < 
oo a.s. Then we have that 


dsI{S*{t)) = 


E* 


ft dsF 2 {u,S*{u),a{u))du 


E* 


d,BSit,S*{t),IiS*it))) 

ff {Fi{u,S*{u),d{u)) + dsFsiu, S*[u),d{u))) du 


2Sd,BS{t,S*{t),I{S*{t))) 
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where 


E, 


= 

2 


+ 


Fi(u, S{u), u(u))du 


e *''G{u,S{u),a{u)) [L'^{u,S{u),a{u)) — a'^{u)) du 

L{u, S{u),a{u))H{u, S{u),(t{u)) (J a^{r)dr \ d[W, M] (u) 


E, 


= -Et 
2 


F 2 {u, S{u), u(u))du 


e *^L'^{u,S{u),a{u)) (j DYcr^{r)drjd[W,M]{u) 


and 


Et 

= -Et 
2 


F3 (u, S{u), a(u))du 


e *'^G{u, S{u),cr{u)) (^L‘^{u, S{u),a{u)) — du 


Proof. See [3] or the previous proof. 


□ 


Remark 8.3. Note that this is generalisation of the formula proved in In 
that case, Fi = F 2 and F 3 = 0 . 


9 Examples 

Now, we give some applications of the decomposition formula for well-known 
models in Finance. 

9.1 Heston Model. 

We consider that the stock prices follows the Heston Model ([T]). Using The¬ 
orem [3H] or Theorem 15.21 we have 


Vit) 


BS{t,X{t),v{t)) 

rT 


+ ^Et 


E* 


e H{u, X{u),v{u)) e ^^drja^{u)vdu 
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Using Theorem 1 7. 11 we have that 


V{t) = BS{t,St,a{t)) 


2 


e H{u, S{u),a{u)) D^a‘^{r)drj a{u)du 

where a^{r) = va{u) exp /J ^dW (s) + 


-k- 


9.2 SABR Model 

We consider that the stock prices follows the SABR model ([2]). Using The¬ 
orem [3T] or Theorem 15.21 we have 


Vit) = BSit,Sit),vit)) 




^-r(u- 


S{u),v{u))a^{u) du 


rE, 


-Et 

2 


A:(m, S{u),v{u))d [M, M] 


(u) 


e *'^H{u,S{u),v{u))a{u)d[W,M]{u) 


where 


and 


d[M,M] =Aa'^a^{t) e“ dt 

d [M, W] = 2aa‘^{t) dt. 


Using Theorem 1 7. 11 we have that 
V{t) = Et[BS{t,Sit),d{t))] 


-Et 

2 




S{u), a{u))a'^{u) — 1^ du 


+ fE. 


e H{u, S{u),a{u)) a'^{r)drj a(u)dh 

where D^a^{r) = 2aCT^(M)l[o,r] (w). 


10 Conclusion 

In this paper, we notice that the idea used in P] can be used for a generic 
Stochastic Differential Equation (SDE). There is no need to specify the volatility 
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process, only existence and uniqueness of the solution of the SDE are needed, 
allowing much more flexibility in the decomposition formula. We see the effect 
on assuming that the stock price follows a exponential process and how a new 
term arises in a general framework. Also, we have computed the decomposition 
using three different method: ltd formula, functional ltd calculus and Malliavin 
calculus. In the case of call options, the idea used in [5] is equivalent to the use 
of the functional ltd formula developed in [3 El 0 [S], but without the need of 
the theory behind the functional ltd calculus. Both formulas can be applied to 
Dirichlet process, in particular, to the fractional Brownian motion with Hurst 
parameter equal or bigger than Furthermore, we realize that the Feynman- 
Kac formula has a key role into the decompossition process. 
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